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ABSTRACT 

Several  integration  schemes  are  applied  to  the  differential  method  for  a sinusoidal- 
profiled  surface-relief  grating  made  of  an  anisotropic  and  conducting  material.  The  nu- 
merical results  show  the  importance  of  the  numerical  stability  and  the  advantage  of  the 
implicit  integration  schemes. 

DIFFERENTIAL  THEORY  OF  GRATINGS 

We  investigate  the  diffraction  problem  on  a surface-relief  grating  ruled  on  an  aniso- 
tropic and  homogeneous  substrate  schematically  shown  in  Fig.l.  The  grating  grooves 
are  parallel  to  the  z -axis  and  the  equation  of  the  grating  surface  is  y = p(x)  where 
p(x ) is  a known  periodic  function  with  the  period  d and  the  depth  h . The  region 
y > p(x)  is  filled  with  a homogeneous  and  isotropic  material  described  by  the  relative 
permittivity  e,  and  the  relative  permeability  p, , and  the  homogeneous  and  anisotropic 
material  that  fills  the  region  y < p(x)  is  described  by  the  relative  permittivity  matrix 

e2  and  the  relative  permeability  matrix  p2  . We  consider  only  time  harmonic  fields  as- 
suming a time-dependence  in  exp(-/co  t) , and  deal  with  the  plane  incident  wave  propa- 
gating in  the  direction  of  polar  angle  0 and  azimuth  angle  <j> . 

The  differential  theory  [1,2]  is  one  of  the  most  commonly  used  approaches  in  the  analy- 
ses of  such  gratings.  Thanks  to  the  periodic  structure,  the  electromagnetic  field  compo 


Figure  1:  Geometry  of  the  surface-relief  grating  under  consideration. 
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nents  can  be  approximately  expanded  in  the  truncated  generalized  Fourier  series  [2];  for 
example  the  x -component  of  E field  can  be  written  as 

Ex(x.y,z)=  Exn  ( v ) exp  | /'  (a  „ x + y - ) ] (1) 


n~  N 


with 


a, 


k]  sin0  cos<(>  + n~~j~ 


y = kt  sin0  sin<J> 


(2) 


where  N is  the  truncation  order.  k]  is  the  wavenumber  in  the  region  y > p(x) , and 
Exll(y ) are  the  n th-order  generalized  Fourier  coefficients  which  are  functions  of  y 

only.  Replacing  all  the  periodic  and  the  pseudo-periodic  functions  by  their  Fourier  se- 
ries and  using  the  Fourier  factorization  rules  [3],  Maxwell's  curl  equations  are  trans- 
formed into  a coupled  differential  equation  set  in  the  form  of 
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where,  for  example,  ex(y)  denotes  a (2/V  + l)xl  column  matrix  generated  by  the  Fou- 
rier coefficients  of  Ex  and  M(y)  is  the  coupling  coefficient  matrix.  Then,  the  scatter- 
ing problem  of  grating  is  reduced  to  an  integration  problem  of  the  coupled  differential 
equation  set  with  boundary  conditions  at  the  top  and  the  bottom  of  the  groove  region. 


NUMERICAL  RESULTS  OF  VARIOUS  INTEGRATION  SCHEMES 

One  method  for  integrating  the  coupled  differential  equation  set  (3)  is  the  rigorous  cou- 
pled-wave method,  which  introduces  the  staircase  approximation  to  describe  arbitrary 
profiled  gratings.  The  real  profile  in  each  step  is  replaced  by  the  structure  uniform  in  the 
y -direction,  and  then  the  boundary-value  problem  can  easily  be  turned  into  an  eigen- 
value problem  because  of  the  absence  of  the  y -dependence.  Another  approach  is  based 
on  the  shooting  method,  which  can  transform  the  boundary-value  problem  into  the  ini- 
tial-value one.  The  initial-value  problem  can  be  solved  by  usual  numerical  integration 
schemes.  In  the  narrow  sense,  this  approach  is  called  the  differential  method.  In  the 
method,  the  Runge-Kutta  or  the  predictor-corrector  Adams  schemes  are  suggested  for 
integration  [2,4], 

Here,  several  numerical  integration  schemes  are  applied  to  the  differential  method  for  a 
sinusoidal  grating  made  of  conducting  material  and  the  convergences  of  the  TM  dif- 
fraction efficiencies  in  -lst-order  with  respect  to  the  number  of  the  integration  steps  are 
compared  in  Fig. 2.  The  grating  parameters  are  chosen  as  follows:  X0  = 0.6328 p/72, 

0 = 30°  , <j)  = 20° , d = 0.6  p/77 , h = 0.5  p/7? , p(x)  = (h/2)[\  + cos(2roc  / d)\ , e,  = p,  = 1 , 

e2..v.v  =£2.w  ~£2~  - -8.19  + /16.38,  e2v;  = — £ 2 -v  =-0.495-/0.106,  e2n.  =£,  n.  = 

e2,,~  =£2,.-i  = 0’  Ei  - I ■ and  TM  (Hv  =0)  polarized  incident  plane  wave.  To  avoid 
the  numerical  difficulty  for  deep  gratings,  we  used  the  scattering  matrix  propagation  al 
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(a)  A = 10 


(b)  N = 25 


Figure  2:  Convergences  of  diffraction  efficiencies  computed  using  various  integration 
schemes  (DM:  differential  method,  RCWM:  rigorous  coupled-wave  method,  ERK: 
classical  fourth-order  Runge-Kutta  scheme,  EAS:  predictor-corrector  Adams  scheme, 
IMS  implicit  midpoint  scheme,  IRK:  implicit  fourth-order  Runge-Kutta  scheme,  IAS: 
implicit  Adams-Moulton  scheme). 


gorithm  [5],  The  rigorous  coupled- wave  method  shows  stable  convergence  but  does  not 
give  reliable  solutions  for  such  a conductive  grating  as  reported  in  Ref.  [4,6].  The  ex- 
plicit schemes  (ERK,  EAS)  show  serious  numerical  instabilities  and  require  large  num- 
ber of  integration  steps  for  reliable  solutions.  On  the  other  hand,  the  implicit  schemes 
(IMS,  IRK,  IAS)  are  numerically  stable  even  when  the  number  of  integration  steps  is 
small.  Consequently,  the  use  of  implicit  schemes  reduces  greatly  the  computation  time 
and  is  highly  suggested. 
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